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Abstract 

We discuss first order transitions for systems in the Ising universality class. The 
critical long distance physics near the endpoint of the critical line is explicitly con- 
nected to microscopic properties of a given system. Information about the short 
distance physics can therefore be extracted from the precise location of the end- 
point and non-universal amplitudes. Our method is based on non-perturbative flow 
equations and yields directly the universal features of the equation of state, without 
additional theoretical assumptions of scaling or resummations of perturbative series. 
The universal results compare well with other methods. 
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1 Introduction 



Many phase transitions are described near the critical temperature by a one- component 
scalar field theory without internal symmetries. A typical example is the water-vapour 
transition where the field <p(x) corresponds to the average density field n(x). At normal 
pressure one observes a first order transition corresponding to a jump in ip from high 
(water) to low (vapour) values as the temperature T is increased. With increasing pressure 
the first order transition line ends at some critical pressure in an endpoint. For p > 
the phase transition is replaced by an analytical crossover. 

This behaviour is common to many systems and characterizes the universality class 
of the Ising model. As another example from particle physics, the high temperature 
electroweak phase transition in the early universe is described by this universality class 
if the mass of the Higgs particle in the standard model is near the endpoint value 
Mjj* ~ 72GeV0. An Ising type endpoint should also exist if the high temperature or 
high density chiral phase transition in QCD or the gas-liquid transition for nuclear matter 
are of first order in some region of parameter space. 

Very often the location of the endpoint - e.g. the critical T*, and n* for the liquid- 
gas transition - is measured quite precisely. The approach to criticality is governed by 
universal scaling laws with critical exponents. Experimental information is also available 
about the non-universal amplitudes appearing in this scaling behaviour. These non- 
universal critical properties are specific for a given system, and the question arises how 
they can be used to gain precise information about the underlying microscopic physics. 
This problem clearly involves the difficult task of an explicit connection between the short 
distance physics and the collective behaviour leading to a very large correlation length. 

So far renormalization group methods 0, |3|, |J have established the structure of this 
relation and led to a precise determination of the universal critical properties. A suitable 
implementation of these ideas should allow us to complete the task by mapping details of 
microscopic physics to non-universal critical quantities. In this paper we show that this is 
indeed possible. A demonstration is given for the liquid-gas transition in carbon dioxide. 

A very useful concept for the study of phase transitions is the coarse grained free 
energy or effective average action T k [ip] |§, §[• It is related to the (grand canonical) 
partition function in presence of an effective infrared cutoff R k , 

Z k [j) = J ^Xexp {-£[*] + J d 3 x (j(x)x(x) - ± x ( x )R k (-d?) X (xf) } , (1.1) 
by a Legendre transform: 

V k [<p] = -lnZ k [j} + Jd 3 x(j{x)i P {x)-^<f{x)R k (-oi 2 )<f(x) 
<p(x) = m ^Z k \j]. 

For the liquid-gas transition x( x ) is the microscopic density field, S[x] = H/T the mi- 
croscopic or classical action and j(x) the source, which will be specified below. The 
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momentum dependent infrared cutoff obeys Rk{q 2 ) ~ k 2 for q 2 <C k 2 and should vanish 
rapidly for q 2 3> k 2 . This assures that only fluctuations with momenta q 2 > k 2 are ef- 
fectively included in the functional integral ( |1 . 1|) . Below, we will often not consider the 
most general form of but rather work with a truncation which only includes the most 
general terms containing up to two derivatives, 

d 3 x {u h (<p(x)) + ^Z k (ip(x))d^d^ . (1.3) 

Our final aim will be the computation of the potential U = U^^o an d the wave function 
renormalization Z = Z^o for a vanishing (or very small) infrared cutoff. In absence 
of the infrared cutoff the shape of U(ip) directly determines the equation of state. For 
a homogeneous situation UT corresponds to the free energy density. Indeed, expressing 
U as a function of the density one finds for the liquid-gas system at a given chemical 
potential /i 

d JL = a. (i.4) 

Equivalently, one may also use the more familiar form of the equation of state in terms 
of the pressure p, 

2 9 (U\ p 



dn\n) T ^ ^ 

(Here the additive constant in U is fixed such that U(n— 0) = 0) . The wave function renor- 
malization Z(ip) contains the additional information needed for the two point correlation 
function at large distance for arbitrary pressure. 

A direct computation of U((p) and Z{tp) becomes, however, a difficult problem for 
all situations where the correlation length becomes large and collective phenomena are 
therefore important. For the water- vapour transition this happens in the vicinity of the 
endpoint, i.e. for a pressure and temperature near p* and T*. In this case the free 
energy can best be calculated by a stepwise procedure where k is lowered consecutively. 
One starts at some large k where only short distance fluctuations are included and a 
relatively simple microphysical computation is reliable. As k is lowered, the effects of 
fluctuations with longer wavelengths are added. We will describe the dependence of the 
average action on k by a non-perturbative flow equation which can be derived from an 
exact renormalization group equation. 

In this way the computation of thermodynamic potentials, correlation length etc. is 
done in two steps: The first is the computation of a short distance free energy T^. This 
does not involve large length scales and can be done by a variety of expansion methods or 
numerical simulations. This step is not the main emphasis of the present paper and we 
will use a relatively crude approximation for the gas-liquid transition. The second step 
is more difficult and will be addressed here. It involves the relation between Ta and r , 
and has to account for possible complicated collective long distance fluctuations. 

For a large infrared cutoff k = A one may compute T\ perturbatively. For example, 
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the lowest order in a virial expansion for the liquid-gas system yields 

(1.6) 



3 MT 

U A (n) = -n ( 1 + Ing + - In ^—^ 



< n \ 6i(A) 2 

+nln {(l-b (A)nW)-— n +CA - 

Here A -1 should be of the order of a typical range of intermolecular interactions, M and g 
are the mass and the number of degrees of freedom of a molecule and bo, b\ parameterize 
the virial coefficient B 2 (T) = b — b\/T . [] (The (mass) density p is related to the particle 
density n by p = Mn.) We emphasize that the convergence of a virial expansion is 
expected to improve considerably in presence of an infrared cutoff A which suppresses the 
long distance fluctuations. 

The field <p(x) is related to the (space-dependent) particle density n(x) by 

<p(x) = KA(n(x)—n) (1.7) 

with h some suitable fixed reference density. We approximate the wave function renor- 
malization Za by a constant. It can be inferred from the correlation length £, evaluated 
at some reference density n and temperature T away from the critical region, through 

dip 2 



r 2 = y,:— 



For a suitable scaling factor 

k a = (- i ; i r ,"z: -~) t <i-9) 




one has Z A — 1. 

We observe that the terms linear in n in eq. ( |1.6|) play only a role for the relation 
between n and p. It is instructive to subtract from Ua the linear piece in ip and to expand 
in powers of ip: 

ua(p) = ^v+^y + ^y+... (i.io) 



with 



mi = K7 2 \^ + 



1 , b (2-b n) 



A \n (l-b n) 2 '~ 

7A = K " (l-b n hY ~v) (L11) 



-r,->/ 68(4-6ofi) | 1' 

1 — Onn 4 ra . 




1 The Van der Waals coefficients bo, b\ for real gases can be found in the literature. These values are 
valid for small densities. They also correspond to k — rather than to k = A. Fluctuation effects lead 
to slightly different values for bi(A) and bi(k = 0) even away from the critical line. We find that these 
differences are small for n<^.n*. Similarly a constant ca should be added to Ua such that t/o(0)=0. 
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For a convenient choice n = — , T = — — one has 7a = and 

36 11 60 

Typical values for carbon dioxide at the endpoint are m A /A 2 = —0.31, Aa/A = 6.63 for 
A -1 = 5 • 10 -10 m, £ — 0.6 A -1 . In the limit (|1.10|) one obtains a tp -model. Our explicit 
calculations for carbon dioxide will be performed, however, for the microscopic free energy 



1761) . The linear piece in the potential can be absorbed in the source term such that the 



equation of state reads Q 



* =j , j = K ?(£ + 1 + * g + ^™). (1 , 3) 



We emphasize that a polynomial microscopic potential ( |1.10 ) with equation of state 



dU/dp = j is a good approximation for a large variety of different systems. For the 
example of magnets ip corresponds to the magnetization and jT to the external magnetic 
field. For 7a = and Aa — * 00, with finite negative m\/\^, this is the ^-symmetric Ising 
model. 

For values of p for which the mass term m 2 {p) = is much larger than A 2 the 

microscopic approximation to Tk remains approximately valid also for k — > 0, i.e. U((p) ~ 
U\((p). The contribution of the long wavelength fluctuations is suppressed by the small 
correlation length or large mass. In the range where m 2 (ip) <C A 2 , however, long distance 
fluctuations become important and perturbation theory looses its validity. Our aim is a 
description of this region by means of non-perturbative renormalization group techniques. 
Beyond the computation of universal critical exponents and amplitude ratios we want to 
establish an explicit connection between the universal critical equation of state and the 
microscopic free energy r A . 

In fig. [l| we plot our results for the equation of state near the endpoint of the critical 
line for carbon dioxide. For the microscopic scale we have chosen A -1 = 0.5 nm. For 
£ = 0.6A -1 , 60(A) = 34 av?mol~ x , 61(A) = 3.11 • 10 6 bar cm 6 mol~ 2 we find the location 
of the endpoint at T* = 307.4 K, = 77.6 bar, p* = 0.442 gcm~ 3 . This compares 
well with the experimental values T* = 304.15 K, = 73.8 bar, p* = 0.468 gcm~ 3 . 
Comparing with literature values 6i(0)w for low density this yields bo(A)/b (0)id = 0.8, 
bi(A)/bi(0)id = 0.86. We conclude that the microscopic free energy can be approximated 
reasonably well by a van der Waals form even for high densities near n*. The coefficients 
of the virial expansion are shifted compared to this low density values by 15-20 per cent. 
The comparison between the "microscopic equation of state" (dashed lines) and the true 
equation of state (solid lines) in the plot clearly demonstrates the importance of the 
fluctuations in the critical region. Away from the critical region the fluctuation effects 
are less significant and could be computed perturbatively. 



2 Note that the source term is independent of k. The linear piece in the potential can therefore easily 
be added to Uk-*o once all fluctuation effects are included. 
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Figure 1: p-/>isotherms of C0 2 for T = 295K, T = 300K, T = 303K, T = 305K, T = 307K 
and T — 315^. The dashed lines represent the virial expansion (at the scale k = A) in 2nd order with 
bo = 34 cm 3 mol^ 1 , b\ = 3.11 • 10 6 bar cm 6 mol~ 2 . The solid lines are the results at the scale k = 

(a = 0.6). 

More details about the equation of state in the critical region as well as an explicit 
expression for U(ri) will be presented in the following sections. In sect. |2| we briefly review 
our method and present the relevant non-perturbative flow equations. Sect. ^| is devoted 
to the critical behaviour of the Ising model. We compute the critical equation of state 
near the endpoint, critical indices, amplitudes and couplings. Our results compare well 
with sophisticated resummation methods of high order renormalization group improved 
perturbation theory, the e-expansion and high temperature series as well as with Monte- 
Carlo simulations and experiment. In this section we also give the universal behaviour of 
the (^-dependent wave function renormalization. This is needed for the dependence of the 
correlation length on density (or a magnetic field in magnets) and usually not available 
by other methods. In sect. f| we turn to the equation of state for first order transitions. 
In the present case of the Ising universality class the universal features can be related to 
the second order phase transition in the Ising model. Such a simple connection is not 
available for more general systems. 

Our flow equation can deal with first order transitions in an arbitrary context. We 
have computed here within a <£> 4 -model with cubic term and for the microscopic free energy 
( |1.6|) . The perfect agreement with the Ising model predictions demonstrates the power 
of this method. The example of carbon dioxide shows the applicability for an arbitrary 
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form of the short- distance potential. 



2 Flow equation for the effective potential 

The dependence of the effective average action on the coarse graining scale k is described 
by an exact functional differential evolution equation || 

It has the form of a renormalization group improved one-loop equation in ((^-dimensional) 
momentum space. The diagonal elements of the effective average propagator, i.e. 

(if } [<p] + Rk)' 1 (q', q) = (T ( k 2) + R k )-\q 2 )(2n) d 6(q' - q) + off-diagonal terms, 



involve the matrix of second functional derivatives of with respect to the Fourier modes 
<p(q), i.e. 5 2 T k [cp)/S<p(-q')5cp(q) = r®[<p]tf,q). 

For a suitable choice of Rk the momentum integral is both infrared and ultraviolet finite 
so that no additional regularization is needed and one can work in arbitrary dimension 
d. This flow equation is formally equivalent to other versions of exact renormalization 
group equations ||. Using instead of the most general form of T the truncated ansatz 
( |1.3| ) one extracts the (approximate) evolution equation for the potential by evaluating 
the flow equation for a constant field <p(q) = (p(2ir) d 5(q), 

T-U k (ip) = v d / dxx* ' 2.2 

at Jo Z k ((p)x + Rk{x) + d A Uk{ip)/d(p 2 

Here t = ln(fc/A), x = q 2 , v^ 1 = 2 d+1 7i2T (d/2), and we choose a smooth infrared cutoff 

R k f x ) = Z ^ kX . (2.3) 

exp(x/A; 2 ) — 1 

The flow equation for Z).(ip) can be obtained by adding to tp a small plane wave with 
long wavelength (see below). We note that the ansatz (|1.3|) can be regarded as the second 
order in a systematic derivative expansion for the solution of eq. (|2.1|). 0, |3j 

Our aim is a numerical solution of the partial differential equation ( |2.2|) with given 
initial conditions at the scale k ~ A. For this purpose it is convenient to introduce a 
dimensionless renormalized field 

<p = k*?Zi£<p (2.4) 
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with wave function renormalization Zo jk = Z k (ipo(k)) taken at the global potential mini- 
mum ifo(k). We also use 

u k (0) = k~ d U k {v) 
z k {(p) = Zo£Z k (<p) 

and denote by u', z' the derivatives with respect to (p. This yields the scaling form of the 
flow equation for u' k : 

J^U£> = ~l( d + 2 - Vo) ■ + \(d-2 + r] )^ ■ ul{{p) ^ 

- 2v d z' k (<p) ■ l d+2 «(£); t?o, h($)) ~ 2 Vdu'k(<f) ■ if «(£); r? , z k (<p)) . 



Similarly, the evolution of z k is described in the truncation (|1.3| ) by 

d _ ..1 

Ofki?) = Vo • z k {(p) +-(d-2 + r}o)<p • z k ((p) 

- - d Vd ' u 'k^) 2 ' m to vo, h(<p)) 
g 

- 2 Vd ' u * ■ m to 2 «(<?); vo, h{0)) ^ ^ 

- -^a-z'^f ■ m to 4 K(£>; Vo, z k {<f>)) ~ 1v d ■ ~z' k \Cp) ■ if m, h{<f)) 
+ ^ d ■ z' k (<p)u'£{(p) ■ l d 2 «((^); r? , z k ((p)) 

+ ~(1 + 2d)v d ■ ~z' k ^f ■ if 2 Vo , ~z k {(p)) 

(see ref. [[J for the derivation of an analogous equation in a different language). Here the 
mass threshold functions 



1 r°° rl 

l d n (u";r] ,z) = - k 2n ~ d -Zl k - dxxz-% 
I ' Jo 



m to{ u "',Vo,z) 



I /"CO A f I 

(P(z) + Z , fe fcV')' 

r 00 , f P 2 (x) 



ljfc 2(„-l)-d . . r dxx iQ 



(2.f 



2 u < fe y i(p(x) + z 0ifc fcvo n . 



(with P(x) = zZ 0tk x + R k (x), P = 4£ and <9 t acting only on P^) describe the suppression 
of fluctuations once Wfc(</5) 3> 1. A discussion of their properties can be found in |7|, [J. 
The anomalous dimension 

Vo, k = - Jt lnZ , k = -Z^-Z^-Z^.— ^.— (2.9) 

is determined by the condition dz(<po)/dt = 0. It appears linearly in the threshold func- 
tions due to d t acting on Z 0jk in the definition ( |2.3| ) of R k . 
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For a computation of 770 we need the evolution of the potential minimum <po(k), which 
follows from the condition -^{dUk/dip(ipo(k)) = 0, namely 

One infers an implicit equation for the anomalous dimension 770,/cj 

4~ _ 8 

Vo = -^d • u'k(Vo) 2 ■ ^4,o( w fc(^o); Vo, 1) + -^d • ul'{<p )z' k (<po) ■ m^+ 2 «(^o); ^0, 1) 

+- d v d ■ 4(^ ) 2 • m^+V^o); *7o, 1) + ■ 4'(^o) • lt«{<fio); Vo, 1) 

-4^ • 4(0oX(^o) • ^2«(^o); 1) (2.11) 

~(l + 2rf)t; d • 4(y5 ) 2 • l d 2 +2 K(Vo); Vo, 1) 

-2«-4H • {z^o)li +2 K^o);vo, i) + <(^o)ifK(^o);%, i)} , 

that can be solved by separating the threshold functions in ^-dependent and ^-indepen- 
dent parts (c.f. eq. ( |2.8|) ). Since ?7o will turn out to be only a few percent, the neglect 
of contributions from higher derivative terms not contained in (|1.3| ) induces a substantial 
relative error for r)o, despite the good convergence of the derivative expansion. We believe 
that the missing higher derivative contributions to 770 constitute the main uncertainty in 
our results. 

For given initial conditions U\((p), Z\(ip) the system of partial differential equations 
(|2.6|) , (|277|) , ( j2. 10| ) , ( |2. 1 ip can be solved numerically. A description of the algorithm used 
for the present work can be found in 0. 



3 Critical equation of state for the Ising model 

In order to make the discussion transparent we present here first results for polynomial 
initial conditions Ql.10 ) with z\ (<p) = 1. The term linear in (p is considered as a source 



j. The special value 7a = realizes the ^-symmetric Ising model. We start with the 
results for the universal critical behaviour for this case. For this particular purpose we 
hold Aa fixed and measure the deviation from the critical temperature by 

8m\ = m\-m\ it = S{T -T c ). (3.1) 
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For the liquid-gas system one has S = 2bi/(K A T 2 ). As predicted by general scaling ar- 
guments we find that for small \Sm\\/A 2 the susceptibility x — m' 1 and the inverse cor- 
relation length or renormalized mass m R = £ _1 in the high temperature phase (5m\ > 0) 
obey 

m 2 = {C + Y l ■ {5m\)\ (3.2) 
m R = (tr'-iSmlr, (3.3) 

whereas for T <T r one has 



m 



2 (C-y 1 ■ (Smlr, (3.4) 



mR = (rr'-iSmlY. (3.5) 
Here the renormalized and unrenormalized squared mass terms are defined by 

d 2 U , , o d 2 U 



m 2 



While for T > T c the minimum of U — Uk-+o is at the symmetric point (fa = 0, the low 
temperature phase is characterized by spontaneous symmetry breaking, with 

<p = B-{-8m\f, (3.7) 

ip% R = 2E-{-5m\y. (3.8) 

The anomalous dimension r\ determines the two point function at the critical temperature 
and equals t\q^ for the scaling solution where d t u = d t z = 0. Our results for the critical 
exponents are compared with those from other methods in table ^. We observe a very 
good agreement for v whereas the relative error for rj is comparatively large as expected. 
Comparison with the lowest order of the derivative expansion (f) shows a convincing 
apparent convergence of this expansion for v. For rj this convergence is hidden by the 
fact that in ]10| a different determination of rj was used. Employing the present definition 



would lead in lowest order of the derivative expansion to a value r\ — 0.11. As expected, 
the convergence of the derivative expansion is faster for the very effective exponential 
cutoff than for the powerlike cutoff (g) which would lead to unwanted properties of the 
momentum integrals in the next order. 

In order to establish the quantitative connection between the short distance param- 
eters m\ and Aa and the universal critical behaviour one needs the amplitudes C ± , 
etc. For A A /A = 5 we find C+ = 1.033, C~ = 0.208, £+ = 0.981, £- = 0.484, 5 = 0.608, 
£' = 0.208. Here and in the following all dimensionful quantities are quoted in units of A. 
The amplitude D is given by dUo/d<p = D ■ ip 5 on the critical isotherme and we obtain 
D— 10.213. In table we present our results for the universal amplitude ratios C + /C~, 
R x = C + DB & -\ R ( = (i+yi»D l '^B. 

The critical exponents and amplitudes only characterize the behaviour of U (<p) in the 
limits if — * ifo and tp — ► oo. Our method allows us to compute U(ip) for arbitrary ip. As 
an example, the quartic coupling X R = |0r(O) = §¥(0), X R = £§-(<Por), p R = \<f% 
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(c) 


0.6300(49) 




1.2400(87) 




(d) 


0.625(1) 






0.025(6) 


(e) 


0.629(3) 






0.027(5) 


(f) 


0.643 


0.336 


1.258 


0.044 


(g) 


0.6181 






0.054 


(h) 


0.6307 


0.3300 


1.2322 


0.0467 


(i) 


0.625(6) 


0.316-0.327 


1.23 - 1.25 





Table 1: Critical exponents of the (d = 3)-Ising model, calculated with various methods. 

(a) Six loop resummed perturbation series at fixed dimension d = 3 ELO]. 

(b) e-expansion in five loop order p|, |l4|], (c) high temperature series ]15| (see also Jl6| , |l7j ) . 
(d),(e) Monte Carlo simulations Jig]- ]2(|. 

(f)-(h): "exact" renormalization group equations. 

(f) effective average action for the 0(iV)-model, N — * 1, with uniform wave function renormalization |l(| 
(see also @). 



(g) scaling solution of equations analogous to (2J3), (2/7) with power like cutoff [|8|. 

(h) effective average action for one-component scalar field theory with field-dependent wave function 
renormalization (present work). 

(i) experimental data for the liquid-vapour system quoted from H 



becomes in the critical region proportional to m^. Our results for the universal couplings 
^B./ m R in the symmetric and A#/m# in the ordered phase can also be found in table [| 
Here tur = |-jr 1^=0 * n the symmetric and rh R = \ ipoR in the ordered phase. 

We should emphasize that the shape of the potential in the low temperature phase 
depends on k in the "inner" region corresponding to \(p\ < <po. This is due to the fluctua- 
tions which are responsible for making the potential convex in the limit k — > | IT| , [12] . [13]] . 
We illustrate this by plotting the potential for different values of k in fig. |2[ 

In the universal critical region the expectation value of (p for very small \Sm\\ and 
sources j (with dU / d(p((fij) = j) can be found from the Widom scaling form of the equation 
of state 

j^U(<p) = M a /(*)-A^ (3.9) 
x = (3.10) 



IV 



Our results for the scaling function f(x) are shown in fig. [3], together with the asymptotic 
behaviour (dashed lines) as dictated by the critical exponents and amplitudes: 

lim /(*) = £> , f(x = -B- 1 /> 8 )=0, 
lim f(x) = (C + )~V , lim f(x) = (Cr^^-a;) 7 (3.11) 
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c + /c- 




R x 










(a) 


4.79 ±0.10 




1.669 ±0.018 






7.88 




(b) 


4.72 ±0.14 




1.648 ±0.031 






9.33 




(c) 


4.95 ±0.015 


1.96 ±0.01 


1.75 






7.9-8.15 




(d) 


4.75 ±0.03 


1.95 ±0.02 








7.76 


5.27 


(f) 


4.29 


1.86 


1.61 


0.865 


0.168 


9.69 


5.55 


(h) 


4.966 


2.027 


1.647 


0.903 


0.204 


8.11 


4.96 


(i) 


4.8-5.2 




1.69 ±0.14 











Table 2: Universal amplitude ratios and couplings of the (e? = 3)-Ising model, 
(a) perturbation theory at fixed dimension d— 3 [JL4J, (b) e-expansion |i~4| . 

(c) high temperature series. Amplitude ratios from B^A/j/mij from ]15|, |l6| , |l7| ]. 

(d) Monte Carlo simulations. Amplitude ratios from po[ , Ar/wr from ]19|, Xr/ttir from Q. 

(f) effective average action for the 0(7V)-model, N — > 1, with uniform wave function renormalization jig ]. 

(h) present work with field-dependent wave function renormalization. 

(i) experimental data for the liquid- vapour system plfl. 



where 

(C-y 1 = B*- l -W±f(-B-W) , f'(x) = I ■ (sf(x) - ^- 5 ) ■ (3-12) 

We have verified the scaling form of the equation of state explicitly by starting with a 
large variety of initial conditions (e.g. different Aa). 

A different useful parameterization of the critical equation of state can be given in 
terms of nonlinearly rescaled fields <pn, using a (^-dependent wave function renormalization 

Z(cp) = Z k ^ ((p), 

<p R = Z{ V f' 2 V . (3.13) 
Our numerical results can be presented in terms of a fit to the universal function 



Fpit(s) = (a s + aj 3 + a 2 s 5 + a 3 s 7 ) ■ f a (s) + (1 - f a (s)) ■ a 4 s 5 . (3.15) 

The factors f a and (1— f a ) interpolate between a polynomial expansion and the asymptotic 
behaviour for large arguments. We use 



/«(*) = a-V- eXPl <?> (3.16) 
l-exp(-Sr) 



A similar fit can be given for 



= linw^j ) s = ^ = /wy /2 = -_ 1/25j (317) 



m)[ 2 V '"/v 
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Figure 2: Uk(<p) for different scales k. The shape of U is stored in smaller intervals At = —0.02 after 
leaving the scaling solution. This demonstrates the approach to convexity in the "inner region" , while 
the "outer region" becomes fc-independcnt. 



zfu(s) = (b + b lS 2 + b 2 s± + b 3 s e + b 4 s 8 )-fp(s) + (l-fp(s))-b 5 \s\-^. (3.18) 

In the symmetric phase one finds (with r\ = 0.0467) a = 1.012, ao = 1.0084, a\ = 3.1927, 
a 2 = 9.7076, a 3 = 0.5196, a 4 = 10.3962 and (3 = 0.5103, b = 1, b x = 0.3397, b 2 = 
-0.8851, 6 3 = 0.8097, b A = -0.2728, 65 = 1.0717, whereas the fit parameters for the 
phase with spontaneous symmetry breaking are a = 0.709, a = —0.0707, a\ = —2.4603, 
a 2 = 11.8447, a 3 = -1.3757, a 4 = 10.2115 and f3 = 0.486, b = 1.2480, 61 = -1.4303, 
62 = 2.3865, 63 = —1.7726, 64 = 0.4904, b 5 = 0.8676 (our fit parameters are evaluated for 
this phase for (d 2 U/dip'j i )(ipR tmax )/k 2 = —0.99). One observes that the coefficients a 2 and 
a 4 are large and of comparable size. A simple polynomial form F = a s + di§ 3 + a 2 s 5 is 
not too far from the more precise result. We conclude that in terms of the rescaled field 
(fin (|3.13| ) the potential is almost a polynomial (^-potential. 

In fig. ^ we show z as a function of s both for the symmetric and ordered phase. Their 
shape is similar to the scaling solution found in ||. Nevertheless, the form of z for k = 
which expresses directly information about the physical system should not be confounded 
with the scaling solution which depends on the particular infrared cutoff. For the low 
temperature phase one sees the substantial dependence of Zk on the infrared cutoff k for 
small values s < sq. Again this corresponds to the "inner region" between the origin 
(s = 0) and the minimum of the potential (sq — 0.449) where the potential finally becomes 
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Figure 3: Widom scaling function f(x) of the (d = 3)-Ising model. ( The present curve is generated 
for a quartic short distance potential Ua with Aa/A = 5). The dashed lines indicate extrapolations of 
the limiting behaviour as given by the critical exponents. 

convex for k — > 0. 

Knowledge of U and z permits the computation of the (renormalized) propagator for 
low momenta with arbitrary sources j. It is given by 



for zq 2 ^ d 2 U/dip 2 R . Here ipn obeys dU/dipn = Zq ' j. We emphasize that the correlation 
length £(<Pr) = z 1 / 2 (Lpn)(d 2 U/d(p 2 R )~" 1 / 2 at given source j requires information about z. 
For the gas-liquid transition is directly connected to the density dependence of 

the correlation length. For magnets, it expresses the correlation length as a function of 
magnetization. The factor z 1 ' 2 is often omitted in other approaches. 

Critical equations of state for the Ising model have been computed earlier with several 
methods. They are compared with our result for the phase with spontaneous symmetry 
breaking in fig. [5] and for the symmetric phase in fig. For this purpose we use F(s) = 
m R 5 ^ 2 dU /difn with s = in the phase with spontaneous symmetry breaking (note 
s ~ s). The constant cp is adapted such that (5= 1) = 1- in the symmetric phase 
we take instead s = J£ ^fi such that ^j(s = 0) = l. One expects for large s an inaccuracy of 

m R 

our results due to the error in rj. 




(3.19) 
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Figure 4: Universal rescaled wave function renormalization z in the symmetric and ordered phase of 
the (d=3)-Ising model. In the low temperature phase the plots are for different k with -£i^p-(ipR. m ax)~c 

(c = —0.9, —0.95, —0.99). Here <pn. m ax is the location of the local maximum of the potential in the inner 
(non-convex) region. In the graph for the high temperature phase we have inserted Monte Carlo results 
by M. Tsypin (private communication). 





2 

m A,crit 


C+ 


D 


A A = 0.1 


-6.4584 • 10~ 3 


0.1655 


5.3317 


A A = 1 


-5.5945 ■ 10" 2 


0.485 


7.506 


Aa = 5 


-0.22134 


1.033 


10.213 


A A = 20 


-0.63875 


1.848 


16.327 



Table 3: The critical values m\ crit and the non-universal amplitudes C + , D as a function of the quartic 
short distance coupling Aa (all values expressed in units of A). Other non- universal amplitudes can be 
calculated from the universal quantities of table ||. 

In summary of this section we may state that the non-perturbative flow equations in 
second order in a derivative expansion lead to a critical equation of state which is well 
compatible with high order expansions within other methods. In addition, it allows to 
establish an explicit connection between the parameters appearing in the microscopic free 
energy T\ and the universal long distance behaviour. For a quartic polynomial potential 
this involves in addition to the non-universal amplitudes the value of m\ crit . We have 
listed these quantities for different values of Aa in table ^. Finally, the temperature scale 
is established by S = dm\/dT \ T . 
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Figure 5: The critical equation of state in the ordered phase. 

(1) mean field approximation. 

(2) effective average action with uniform wave function renormalization ]io|| . 

(3) Monte Carlo simulation Jl9| . 

(4) resummed e-expansion in 0(e 3 ), five loop perturbative expansion and high temperature series |l 

(5) present work. 



4 Equation of state for first order transitions 

Our method is not restricted to a microscopic potential with discrete ^-symmetry. The 
numerical code works for arbitrary initial potentials. We have investigated the polynomial 
potential (|1.10p with 7a ^ 0. The numerical solution of the flow equations (|2.6|) , (|2.7f) 
shows the expected first order transition (in case of vanishing linear term j). Quite 
generally, the universal critical equation of state for first order transitions will depend on 
two scaling parameters (instead of one for second order transitions) since the jump in the 
order parameter or the mass introduces a new scale. The degree to which universality 
applies depends on the properties of a given model and its parameters (see || [TIJ for a 
detailed discussion). For a <£> 4 -model with cubic term ( 1.10| ) one can relate the equation 



of state to the Ising model by an appropriate mapping. This allows us to compute the 
universal critical equation of state for arbitrary first order phase transitions in the Ising 
universality class from the critical equation of state for the second order phase transition 
in the Ising model. For other universality classes a simple mapping to a second order 
equation of state is not always possible - its existence is particular to the present model. 
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Figure 6: The critical equation of state in the symmetric phase. 
(1) Monte Carlo simulation [p2| , (2) e-expansion jwj. 

(3) effective average action with uniform wave function renormalization ]To| . 

(4) Monte Carlo simulation |Q, (5), (6) high temperature series [pT[,]l6[. 
(7) present work. 



By a variable shift 



a 



7a 
3A a 



we can bring the short distance potential ( |1.10| ) into the form 

A 



Z o 



with 



/4 



IK 

3A 



27Ai 



m A - 



6A a 



(4.1) 



(4.2) 



(4.3) 



We can now solve the flow equations in terms of a and reexpress the result in terms of ip 
by eq. ([4.1|) at the end. The exact flow equation (|2.1| ) does not involve the linear term 
~ J 7 cr on the right hand side (also the constant c\ is irrelevant). Therefore the effective 
potential {k = 0) is given by 



U = U Za (a)-J~<r + c A = U Z2 (ip + ^)-(ip + ^)J y + c A 



7a 



3A/ 



3A/ 



(4.4) 
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J + Jy, (4-5) 



where U Z2 is the effective potential of the Ising type model with quartic coupling Aa and 
mass term /i A . The equation of state dU/dtp = j or, equivalently 

dU Z2 

is therefore known explicitly for arbitrary m A , 7a and Aa (cf. eq. ( |3.9| ) for the universal 
part). This leads immediately to the following conclusions: 

i) First order transitions require that the combination U ((f) — jip has two degenerate 
minima. This happens for J 7 + j = and fi\ < /z A crit or 

7a ( 2 7a \ U c\ 

m A - — (4.6) 



3A A V 9A A/ 

2 

< < lA,crit + 7T-- ( 4 - 7 ) 



6A A 

Here /i A crit is the critical mass term of the Ising model, 
ii) The boundary of this region for 

6A A 

is a line of second order phase transitions with vanishing renormalized mass or 
infinite correlation length. 

For j = (e.g. magnets with polynomial potential in absence of external fields) the 
equations (|4.6| ), (|4.8| ) have the solutions 

7a ; i = , 7 A;2,3 = ±(-18A A ^l c „ t ) 1/2 . (4.9) 

The second order phase transition for 7a 7^ can be described by Ising models for shifted 
fields a. For a given model, the way how a phase transition line is crossed as the tem- 
perature is varied follows from the temperature dependence of j, m\, 7a and Aa- For the 
gas-liquid transition both j and m\ depend on T. 

In the vicinity of the boundary of the region of first order transitions the long range 
fluctuations play a dominant role and one expects universal critical behaviour. The de- 
tailed microscopic physics is only reflected in two non-universal amplitudes. One reflects 
the relation between the renormalized and unrenormalized fields as given by Z$. The 
other is connected to the renormalization factor for the mass term. Expressed in terms of 
renormalized fields and mass the potential U looses all memory about the microphysics. 

The critical equation of state of the non-symmetric model (7A 7^ 0) follows from the 
Ising model (4.5). With |^ = \ip\ s f(x), the scaling form of the equation of state 
j = ?F ^ or ^ ne m °d e l with cubic coupling can be written as 



1 . 7a , s ,, v / 7a 2 7a \ 7a r 2 (a 1n x 

3 = l(p+ VTj f(x) - 3AT^ + 54AT -3A7^ A ' (410) 
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where ./• = U/R and 5[i\ = m\ — — ^\ crit - One may choose 

„2 



6A A 



V = ^i^+l^ + ^J^+SAl 1 (411) 

as the second scaling variable. For small symmetry breaking cubic coupling 7a one notes 
y ~ 7a. The scaling form of the equation of state for the non-symmetric model reads 

3 = + {/(*)- ?/}. (4.12) 

This universal form of the equation of state is relevant for a large class of microscopic free 
energies, far beyond the special polynomial form used for its derivation. 

It is often useful to express the universal equation of state in terms of renormalized 
fields and masses. We use the variables 

<Pr rn R 
s = — ; v = — gj, 4.13 
Vor m% v ' 

where m# = (f^r l^ 0J? J is the renormalized mass at the minimum <^ i? of U(^pr) 

whereas m^ 2 is the renormalized mass at the minimum of the corresponding ^-symmetric 
effective potential obtained for vanishing cubic coupling 7a = 0. Then the critical tem- 
perature corresponds to v — 1. In this parameterization the universal properties of the 
equation of state for the Ising type first order transition can be compared with transitions 
in other models - e.g. matrix models M - where no simple mapping to a second order 
phase transition exists. 

A convenient universal function G(s, v) for weak first order transitions can be defined 

as 

G(»,«):=^. (4.14) 

We plot G(s, v) in fig. [7] as a function of s for different values of v. For the present 
model all information necessary for a universal description of first order phase transitions 
is already contained in eqs. ( f4.4|) or ([4.5| ). The function G(s,v) can serve, however, for a 
comparison with other models. At this place we mention that we have actually computed 
the potential U both by solving the flow equations with initial values where 7a ^ 
and by a shift from the Ising model results. We found good agreement between the two 
approaches. 



5 Discussion 

In conclusion, we have employed non-perturbative flow equations in order to compute 
explicitly the equation of state. We have first studied models where the microscopic free 
energy can be approximated by a polynomial approximation with terms up to quartic 
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Figure 7: G(s, v) for v = 1, v = 1.037, u = 1.072, v = 1.137 and v = 1.195. 



order. This covers second order as well as first order transitions, both for the universal and 
non-universal features. The same method can be used away from the critical hypersurface, 
allowing therefore for an explicit connection between critical and non-critical observations. 

The ability of the method to deal also with a microscopic free energy which is not 
of a polynomial form is demonstrated by a particular example, namely the equation of 
state for carbon dioxide. In the vicinity of the endpoint of the critical line we can give an 
explicit formula for the free energy density U(n,T)T. Using the fits (3.15), ((3.18), for 
F(s) and z(s) one finds 3 



3 Note that is somewhat different form n and therefore <pr is defined slightly different from eq. 
(1.7). This variable shift (similar to (4.1) reflects the fact that eq. (1.6) contains higher than quartic 
interactions and cannot be reduced to a </? 4 -potential even for 7a = 0. 



U(n,T) = U Z2 (Mn,T), m R (T)) + J(T)(n - n*) - K(T) 




(5.2) 



with fij ~ Oj and 
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<f R (n, T) = H± 
m R (T) = £ T ± 



rp rp — r\V 

* ( \ 

[n — n* J 



T-T* 



(5.3) 



The two non-universal functions J(T) and K(T) enter in the determination of the chemical 
potential and the critical line. 

In particular, the nonuniversal amplitudes governing the behaviour near the endpoint 
of the critical line can be extracted from the equation of state: In the vicinity of the 
endpoint we find for T = T* 



P>-P* = P*-P< = DlUV-™) (5.4) 




with Dp = 2.8 g~ l cm i 1 where p> > p* and p< < p* refer to the density in the high and 
low density region respectively. At the critical temperature T c < T* and pressure p c < p* 
for a first order transition one finds for the discontinuity in the density between the liquid 
(pi) and gas (p g ) phase 




Ap = p.-P, = BAZ-JZ) =B T (^—^\ (5.5) 



with Bp = 0.85 gcm~ 3 , = 1.5 gcm~ 3 . This relation also defines the slope of the critical 
line near the endpoint. 

There is no apparent limitation for the use of the flow equation for an arbitrary 
microscopic free energy. This includes the case where U\ has several distinct minima and, 
in particular, the interesting case of a tricritical point. At present, the main inaccuracy 
arises from a simplification of the g 2 -dependence of the four point function which reflects 
itself in an error in the anomalous dimension rp The simplification of the momentum 
dependence of the effective propagator in the flow equation plays presumably only a 
secondary role. In summary, the non-perturbative flow equation appears to be a very 
efficient tool for the establishment of an explicit quantitative connection between the 
microphysical interactions and the long-range properties of the free energy. 

Acknowledgment: We would like to thank B. Bergerhoff, J. Berges and M. Tsypin for 
helpful discussions. 
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